Transport properties of graphene quantum dots 
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In this work we present a theoretical study of transport properties of a double crossbar junction 
composed by segments of graphene ribbons with different widths forming a graphene quantum 
dot structure. The systems are described by a single-band tight binding Hamiltonian and the 
Green's function formalism using real space renormalization techniques. We show calculations of 
the local density of states, linear conductance and I-V characteristics. Our results depict a resonant 
behavior of the conductance in the quantum dot structures which can be controlled by changing 
geometrical parameters such as the nanoribbon segments widths and relative distance between them. 
By applying a gate voltage on determined regions of the structure, it is possible to modulate the 
transport response of the systems. We show that negative differential resistance can be obtained for 
low values of gate and bias voltages applied. 



PACS numbers: 61.46.-w, 73.22.-f, 73.63.-b 
I. INTRODUCTION 

In the last few years, graphene-based systems have 
attracted a lot of scientific attention. Graphene is a 
single layer of carbon atoms arranged in a two dimen- 
sional hexagonal lattice. In the literature, it has been 
reported several experimental techniques in order to ob- 
tain this crystal, such as mechanical peeling or epitax- 
ial growtlP^. On the other hand, graphene nanoribbons 
(GNRs) are stripes of graphene which can be obtained by 
different methods like high-resolution lithograph};^, con- 
trolled cutting processes^ or by unzipping multiwalled 
carbon nanotubes^. Different graphene heterostructures 
based on patterned GNRs have been proposed and con- 
structed, such as graphene junction^, graphene flaked, 
graphene antidots superlattice^, and graphene nano- 
constrictions^^. The electronic and transport prop- 
erties of these nanostructures are strongly dependent 
of their geometric confinement, allowing the possibil- 
ity to observe quantum phenomena like quantum inter- 
ference effects, resonant tunneling and localization. In 
this sense, the controlled modification of these quan- 
tum effects by means of external potentials which change 
the electronic confinement, could be used to develop 
new technological applications such as graph ene-ba sed 
composite material^^l^ molecular sensor device^^^^ and 
nanot r ansist or s ^ ^. 

In this work we study the transport proper- 
ties of quantum-dot like structures, formed by seg- 
ments of graphene ribbons with different widths con- 
nected between each other, forming a double crossbar 
junctior^^^. These graphene quantum dots (GQDs) 
could be versatile experimental systems which allow a 




FIG. 1: Schematic view of a GQD structure based on leads 
of width Nl = 9, and a conductor region composed by two 
symmetrical junctions of width Nb — 21 and length Lb — ^ 
separated by a central structure of length Lc = 4 and width 
Nc = 9 



range of operational regimes from conventional single- 
electron detectors to ballistic transport. The systems we 
have considered are conductors formed by two symmet- 
ric crossbar junctions of widths Nb and length L^, and 
a central region that separates the junctions, of width 
Nc and length Lc- Two semi-infinite leads of width 
Nl = Nc are connected to the ends of the central con- 
ductor. A schematic view of the considered system is 
presented in figure [l] We studied the different electronic 
states appearing in the system as a function of the ge- 
ometrical parameters of the GQD structure. We found 
that the GQD local density of states (LDOS) as a func- 
tion of the energy shows the presence of a variety of sharp 
peaks corresponding to localized states and also states 
that contributes to the electronic transmission which are 
manifested as resonances in the linear conductance. By 
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changing the geometrical parameters of the structure, it 
is possible to control the number and position of these 
resonances as a function of the Fermi energy. On the 
other hand a gate voltage applied at selected regions of 
the conductor allows the modulation of their transport 
properties exhibiting a negative differential conductance 
(NDC) at low values of the bias voltage. 

II. MODEL 

All considered systems have been described by using 
a single 7r-band tight binding hamiltonian, taking into 
account nearest neighbor interactions with a hopping 
parameter 70 = 2.75 eV. Besides, we have considered 
hydrogen passivation by setting a different hopping pa- 
rameter for the carbon dimers at the ribbons edges^^, 

ledge = 1.127o . 

The electronic properties of the systems have been cal- 
culated using the surface Green's functions matching for- 
malism (SGFMF)P^I^. In this scheme, we divide the 
heterostructure in three parts, two leads composed by 
semi-infinite pristine GNRs, and the conductor region 
composed by the double GNRs crossbar junctions, as it 
is shown in figure [l] 

In the linear response approach, the electronic conduc- 
tance is calculated by the Landauer formula. In terms of 
the conductor Green's functions, it can be written a^^ 

G=^f (E) = ^Tr [r.G^rnG^] , (1) 

where T{E), is the transmission function of an electron 
crossing the conductor region, T^/ji = ^[I]^,/^ — S^^^^J 
is the coupling between the conductor and the respec- 
tive leads, given in terms of the self-energy of each lead: 
^L/R = Vc,L/RgL/RVL/R,C' Here, Vc,l/r are the cou- 
pling matrix elements and Ql/r is the surface Green's 
function of the corresponding leacP^. The retarded (ad- 
vanced) conductor Green's function are determined b}E3 

G^'^ = [E-Hc-Ef'' -E^'""]-' (2) 

where He is the hamiltonian of the conductor. In order 
to calculate the differential conductance of the systems, 
we determine the I-V characteristics by using the Lan- 
dauer formalism^. At zero temperature, it reads 

/io+V/2 

/(y) = ^ I f{E,V)dE, (3) 

where jhq is the chemical potential of the system in equi- 
librium and T{E,V) is defined by equation ([l]). The 
Green's functions and the coupling terms depend on the 
energy and the bias voltage. We consider a linear volt- 
age drop along the longitudinal direction of the conductor 
and the gate voltage is included in the on-site energy at 



the regions in which this potential is applied. In what fol- 
lows the Fermi energy is taken as the zero energy level, 
all energies are written in terms of the hopping param- 
eter 7o and the conductance is written in units of the 
quantum of conductance Gq = 2e^ /h. 



III. RESULTS AND DISCUSSION 

In figure |2j we display results of the linear conductance 
for a graphene quantum dot structure formed by two 
armchair ribbons leads of width = 5 and a conductor 
region composed by two symmetric crossbar junctions of 
width Nb = 17 and variable lengths Lb (from 1 up to 7) 
. Two relative distances between the junctions: panel (a) 
Lc = 5 and panel (b) Lc = 10 are considered and the 
conductance of a pristine Nl = 5 armchair nanoribbon 
is included as a comparison (light green dotted line). 

In both panels it is possible to observe a series of peaks 
at defined energies in the conductance curves. This res- 
onant behavior of the electronic conductance arises from 
the interference of the electronic wave functions inside 
the structure, which travel forth and back forming sta- 
tionary states in the conductor region (well-like states). 
In order to understand these results, it is convenient to 
define two energy regions, the low-energy range from 
up to 0.7 7o (corresponding to the first quantum of con- 
ductance for the pristine N = 5 armchair ribbon) and 
the high-energy range, from 0.7 to 1.2 70 (corresponding 
to the second step of conductance of the = 5 pristine 
system). 

In the low-energy range, it is clear that the conduc- 
tance peaks correspond to resonant states belonging to 
the central region of the conductor. By increasing the 
relative distance Lq of the central part of the system, 
the number of allowed well-like states also increases and, 
as a cons equence , the conductance curves exhibit more 
resonances^^'^^^. The well-like states remain almost in- 
variant under geometrical modifications of the crossbar 
junctions. However, for certain energy ranges and for 
particular junction lengths, the electronic transmission of 
the systems, exhibits an almost constant value. For in- 
stance, in both panels of figure [2j for the cases of Lb = 1 
and Lb = 4 at the energy range 0.4 to 0.65 70. This effect 
corresponds to a constructive interference between well- 
like states from the central region with states belonging 
to the crossbar junctions regions. The different inter- 
ference effects will be clarified by analyzing the LDOS 
of these systems, which is shown next in this paper. In 
the high-energy region, the conductance curves exhibit a 
complex behavior as a function of the geometrical param- 
eters of the GQD structures. There is not a predictable 
behavior of the conductance as the width and length of 
the crossbar junction are increased. 

It is important to pointed out, from the analysis of 
figure |2j that it is possible to identify some interesting 
effects associated to well known quantum phenomena. 
For instance, in panels (a) and (b) of figure [2j for the 
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FIG. 2: Conductance as a function of the Fermi energy for a 
graphene quantum dot structure composed by two armchair 
ribbons leads of width Nl = 5, a double symmetric crossbar 
junctions of width Nb = 17 and variable length, from Lb — ^ 
up to Lb = 7. The central region has a width Nc = 5 
and two separations (a) Lc = 5 and (b) Lc = 10. Light 
green dotted line corresponds to the conductance of a pristine 
Nl = 5 ribbon. All curves have been shifted 2Go for a better 
visualization. 



cases Lb = 5,6 and 7 at energies around E = O.570, it 
is possible to observe a non- symmetric line shape, which 
corresponds to a convolution of a Fano-lik^^and a Breit- 
Wignei^^ resonances. This kind of line-shape, has been 
observed before in other mesoscopic systems by Orellana 
and co-worker^^. In that reference, a simple model of 
two localized states with the same energy ujc are directly 
coupled between each other by a coupling r and indi- 
rectly coupled throughout a common continuum. The 
corresponding resonances have been adjusted by using 
the following expression: 



T{oj) 




[(1 - a;2)r/2 - (w - a;,)2 + r^]2 + _ _ rx)^ 

(4) 

where 77 is the width of a localized state coupled to the 



FIG. 3: (Color online) Numerical adjustment of a convolution 
of a Fano and a Breit -Wigner line shapes (red dash line) and 
one of the conductance resonance in the system (black solid 
line) with r] = O.OI70, cJc = 0.597o, x=0.7,t = O.O270. 



continuum and x defines the degree of asymmetry of the 
system. 

We realize that a possible interference mechanism oc- 
curring in our considered system can be explained with 
the above model, which helps to get an intuitive under- 
standing of the origin of some conductance lines shape. 
In figure [3] we have plotted a particular conductance res- 
onance and the corresponding fitting given by the model 
represented by equation Q , where it is observed the good 
agreement between both curves. 

In what follows we focus our analysis on the resonant 
behavior exhibited by the conductance curves, analyz- 
ing the different electronic states in the conductor. We 
have performed calculations of the spatial distribution 
of LDOS for certain energies corresponding to different 
states present in the conductor. In the bottom panel of 
figure [4] we show results for the LDOS as a function of the 
Fermi energy, for a GQD structure formed by a double 
crossbar junction of width Nb = 17 and length Lb = 
separated by a central region of width Nq = 5 and length 

We start our analysis focusing on the some sharp states 
present in the curve LDOS vs energy of this figure. We 
have marked the first three sharp states in this LDOS 
plot with the letters (a), (b) and (c) and we have calcu- 
lated the spatial distribution of these states, representing 
by the corresponding contour plots exhibited in the fig- 
ure. These states are completely localized at the crossbar 
junctions, and as we establish in a previous work^, their 
correspond to bound states in the continuum (BICs^P^M^. 
It is not expected that these kind of states play a role in 
the electronic transport of these GQD structures, which 
is shown at the corresponding conductance curves of fig- 
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FIG. 4: LDOS for a GQD formed by a double crossbar junc- 
tion of width Nb= 17 and length Lb— separated by a cen- 
tral region of width Nc= 5 and length Lc=5. Panel (a), 
(b) and (c) correspond to the contour plots of some sharp 
LDOS resonances marked in the bottom plot. As a reference, 
the LDOS of a pristine N=5 armchair ribbon is plotted as a 
dotted green line 

ure[2l 

In what follows we will focus our analysis to those 
states that contribute to the conductance of the systems. 
In figure |5] it is shown the spatial distribution of LDOS 
for a GQD formed by a double crossbar junction of width 
Nb = 17 and length Lb = 3, separated by a central re- 
gion of width Nc = 5 and length Lc = 5. As a reference, 
at the bottom panel we have included a plot with the 
LDOS versus Fermi energy of the GQD system, there we 
have marked with letters (a), (b), (c) and (d) four partic- 
ular energy states. The corresponding contour plots are 
displayed at the upper part of the figure. 

In these plots, it is possible to observe the resonant be- 
havior of these states, which are completely distributed 
along the GQD structure, presenting a maximum of the 
probability density at the center of the system. This con- 
dition favors the alignment of the electronic states of the 
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FIG. 5: LDOS for a GQD structure formed by a double cross- 
bar junction of width Nb= 17 and length Lb= S separated 
by a central region of width Nc= 5 and length Lc=5. Panel 
(a),(b), (c) and (d) correspond to the contour plots of those 
resonant states marked in the LDOS plot displayed at the 
bottom. 



leads with the resonant states in the conductor and con- 
sequently, a unitary transmission at those energy values 
is expected. This behavior is reflected as a series of reso- 
nant peaks in the conductance of the system that could 
be controlled by means of the geometrical parameters of 
the GQD, as it is shown in figure |2j At higher energies 
there is an interplay between localized states in the cross- 
bar junctions and resonant states in the central region of 
the conductor. As it is shown in panel (d) of figure [sj 
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FIG. 6: Conductance as a function of Fermi energy and gate voltage for GQDs composed by leads of width Nl = 5, two 
crossbar junctions of width Nb = 17, a central part of width Nc — 5 and length Lc — 5. Panels (a) and (c) correspond to 
junctions of length Lb = "2, while panels (b) and (d) correspond to junctions of length Lb = S. In the upper panels the gate 
voltage is applied at the crossbar regions, and in the lower panels the gate voltage is applied at the central structure. The black 
segmented lines highlight different slopes discussed in the text 



some states are strongly dependent of the geometry of 
the junctions, therefore for some particular configuration 
it is possible to observe a non-null transmission at these 
energies, while for other cases, there are destructive inter- 
ference mechanisms that suppress completely the trans- 
mission at that energy region. 

We have studied different configurations of GQDs 
structures, varying systematically some geometric pa- 
rameters. We have observed a quite general behavior of 
such resonant conductors with the presence of sharp and 
resonant states. Depending of each particular system, it 
can be observed changes in the number and position in 
energy of these states, as well as in their spatial distri- 
bution. The different intensity of the peaks in the LDOS 
curves depends on the spatial distribution of the states. 
There are states completely extended along the conduc- 
tor (like in panel (b) of figure 5] ) which generate wider 
and less intense peaks. On the other hand, there are 
resonant states which are more concentrated in certain 
regions of the conductor (like as panel (c) and (d) of fig- 
ure [5| which generate sharper and more intense peaks in 
the LDOS. 

Now we focus our analysis on the effects of an applied 
gate voltage on the transport properties of GQD struc- 
tures. Results of the conductance as a function of the 
Fermi energy, for different values of a gate voltage ap- 
plied in selected regions of a GQD are shown in figure 
[6] The systems are composed by leads of width Nl = 5, 
two crossbar junctions of width Nb = 17, a central part 
of width Nc = 5 and length Lc = 5. Panels (a) and (c) 
correspond to junctions of length Lb = 2^ while panels 
(b) and (d) correspond to junctions of length Lb = Fi- 
nally, upper panels correspond to a gate voltage applied 
at the crossbar junctions regions, whereas the lower pan- 



els correspond to a gate voltage applied at the central 
part of the structure. 

In these contour plots of conductance, it is possible to 
observe the behavior of the resonant states of the sys- 
tem with a gate voltage applied at different regions of 
the structure. The lower panels of the figure [6] show the 
case of a gate voltage applied at the central region of the 
considered GQDs. In these plots the linear dependence 
of the conductance resonances as a function of the gate 
voltage is manifested. It can be shown that the electronic 
states of a pristine armchair graphen e ribb on are regu- 
larly spaced in the whole energy rang^^^, therefore, as 
the gate voltage is applied, there will be a high probabil- 
ity that the lead states become aligned with the resonant 
states in the central region of the structure. This behav- 
ior in completely general and independent of the width 
Lb of the crossbar junctions. The linear behavior of the 
conductance peaks could be useful in nanoelectronics de- 
vices, due to the possibility of control the current flow 
through these systems. This argument will be more clear 
with the analysis of the differential conductance, which 
is shown below in this paper. 

The case of a gate voltage applied at the crossbar junc- 
tion regions is shown in the upper panel of figure [6| The 
conductance behavior in more complicated to analyze, 
nevertheless it is still possible to observe a linear de- 
pendence of the resonant states with the gate voltage. 
However, two different slopes can be noticed, for states 
belonging to the crossbar junctions (lower slope) and for 
states belonging to the central region of the conductor 
(higher slope). Besides, the panels exhibit an important 
reduction in the conductance gap, for different values of 
gate voltage. This effect is mainly produced by an en- 
ergy shift of the localized states at the junctions, which 
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induces a less destructive interference with the resonant 
states. It is important to point out that this effect only 
can be observed because the gate potential is applied si- 
multaneously to both crossbar junctions, otherwise, the 
conductance gap would be not noticeably affected by the 
gate potential. Finally, the dark (blue online) regions 
present in figure |6j occur at energy ranges around the 
LDOS singularities of the pristine N=5 AGNR. At these 
energies, it appear the second and the third allowed states 
appear, which interrupt the linear behavior of the con- 
ductance resonances as a function of the gate voltage. 

In order to understand the presence of two different 
slopes in the upper panel of figure |6) we present a simple 
model which keeps the underlying physics of the con- 
sidered system and allow us to explain qualitatively our 
results. 
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FIG. 7: A simple model of two crossbar junctions formed by 
two quantum dots, coupled to a linear chain of sites. 

The scheme showed in figure [7| is a simple represen- 
tation of our conductor. The system is composed by a 
linear chain of three sites, which are connected to two 
semi-infinite leads. We have considered four quantum 
dots connected to the extremes of the chain forming a 
double crossbar junction, at which we have applied sym- 
metrically a gate voltage Vg. This potential will modify 
the on-site energy of the dots by a linear shift of energy 
proportional to the gate voltage amplitude. 

By using Dyson equation, it is possible to calculate the 
Green's function of the central site of the chain labeled 
by 0, which takes the form: 



00 



- So 



(5) 



where u is the energy of the incident electrons, is the 
central on-site energy and the self energy S is given by 
the following expression (see Appendix for a detailed de- 
duction): 



2 _ 

+ r2 L 



Vn 



(6) 



In this model, the self energy of the Green's functions 
of the central region acquires a real part that depends 
on the gate voltage applied to the crossbar junction re- 
gion. As a consequence, two different slopes appear in 



the behavior of the conductance peaks as a function of 
the gate voltage. One of these slopes corresponds to the 
direct evolution of the states belonging to the crossbar 
junctions as a function of the gate voltage (lower slope), 
and the higher slope corresponds to the indirect states 
belonging to the central region. 
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FIG. 8: (a) Current versus bias voltage and (b) Differential 
conductance as a function of bias voltage for a GQD composed 
by two crossbar junctions of length Lb = 2 and width Nb — 
17 separated by a central region of width Nc = 5 and a 
variable length from Lc = 1 up to Lc = 20. All curves have 
been shifted for a better visualization. 

Now we focus our analysis on the I-V characteristics 
and the differential conductance of these resonant GQDs. 
Figure [S] shows results of these transport properties for 
a conductor formed by two crossbar junctions of width 
Nb = 17, length Lb = ^2 and a central region of width 
Nc = 5. In this plots the length of the central struc- 
ture is varied from Lc = 1 up to Lc = 20. In panel (a) 
it is possible to observe that for a very small separation 
between both junctions, the I-V characteristics exhibit 
abrupt slope changes, and oscillations for certain ranges 
of the bias voltage. This behavior is produced by the 
increasing number of resonant states as a result of the 
enlargement of the conductor central region. The ap- 
plied bias voltage allows the continuum alignment of the 
resonant states of the system with the electronic states of 
the leads, leading to variations of the current intensity. 
On the other hand, every I-V curve shows a wide gap 
of zero current until certain bias voltage. This threshold 
value exhibits a linear dependence of the length of the 
central region of the conductor. As the relative distance 
between the junctions is increased, there are more reso- 
nant states available at lower energies because the elec- 
tronic confinement in this region is weaker, therefore, the 
electronic transmission under a bias voltage is possible at 
lower voltage values. 

The abrupt changes of the current as a function of 
the bias voltage are clearly reflected in the differential 
conductance of these systems. In panel (b) of figure [sjit is 
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possible to observe the behavior of the dl/dV — V curves 
as a function of the length of the central region of the 
conductor. As this region is enlarged, the oscillations in 
the differential conductance become more evident. Each 
time the bias voltage aligns the resonant states of the 
conductor with leads states, the current will increase and 
a positive change in the differential conductance occurs. 
However, if the bias voltage is not enough to align the 
states, the current drops and the differential conductance 
becomes negative in a range of voltage. This can be seen 
in the cases Lc = 10, Lc = 15 and Lc = 20 in the figure 
[s] (b). The bias voltage value at which occurs negative 
differential conductance (NDC)^^ depends directly on 
the relative distance between the crossbar junctions. 
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FIG. 9: Differential conductance of a GQD structure com- 
posed by two crossbar junctions of length = 2 and width 
Nb = 17 separated by a central region of width Nc = 5 and 
length Lc — 5. The gate voltage is applied at selected regions 
of the conductor: (a) on the crossbar junction region and (b) 
on the central region. 

Now we focus our analysis in the effects of a gate volt- 
age on the differential conductance of the systems. In 
figure |9] we present results of the differential conductance 



of a GQD composed by two crossbar junctions of length 
Lb = 2 and width Nb = 17 separated by a central region 
of width Nc = 5 and length Lc = 5. The gate poten- 
tial has been applied at selected regions of the structure: 
(a) crossbar junction regions and (b) at the central re- 
gion. In panel (b) of this figure it can be observed a 
periodic modulation of the differential conductance as a 
function of the gate voltage applied at the central region 
of the GQD. This behavior is determined and directly 
related with the linear evolution of the resonant states 
of the conductor, and consequently, with the peaks of 
conductance of the system [panel (c) and (d) of figure 
[6]. For instance, in the configuration considered in this 
figure, it can be noted abrupt changes of the differential 
conductance as a function of the gate voltage for bias 
voltages values around O.870, which indicates an abrupt 
increasing of the current fiowing through the conductor. 
This behavior is very general for other systems studied, 
which suggests possible applications in the development 
of GQDs based electronic devices. 

In the case in which the gate voltage is applied simul- 
taneously at the crossbar junction regions, (panel (a) of 
this figure) it is still possible to note certain regularity in 
the dependence of differential conductance as a function 
of the external potentials, although a periodic modula- 
tion is not clearly observed. However, in this configura- 
tion of applied potentials our results show that the GQD 
structure exhibits NDC at very low values of bias volt- 
age and gate voltages. This behavior can be understood 
by observing panels (a) and (b) of figure |6j where the 
contour plots show areas where the conductance is com- 
pletely suppressed at low energies, for different values of 
gate voltage. We have observed this kind of behavior in 
every configuration considered in this work. 

In relation to the practical limitations of our calcula- 
tions we would like to mention that despite of the size of 
the structures used in this work are below the limit of the 
actual experimentally realizable, our calculations can be 
scaled to structures of bigger sizes. By other hand our 
model does not include disorder or electron-electron in- 
teraction, nevertheless we are convinced that our results 
will be robust under these kind of effects as it is in meso- 
scopic systems. For instance it is known that in quantum 
dots, the resonant tunneling and the Fano effe ct sur vive 
to the effect of the electron-electron interactioil^^MI^ 



IV. SUMMARY 

In this work we have analyzed the transport proper- 
ties of a GQD structure formed by a double crossbar 
junction made of segments of GNRs of different widths. 
We have focused our analysis on the dependence of the 
electronic and transport properties with the geometrical 
parameters of the system looking for the modulation of 
these properties through external potentials applied to 
the structure. Our results depict a resonant behavior of 
the conductance in the quantum dot structures which can 
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be controlled by changing geometrical parameters such as 
nanoribbon widths and relative distance between them. 
We have explained our results in terms of the analysis of 
the different electronic states of the system. The possi- 
bility of modulating the transport response by applying 
a gate voltage on determined regions of the structure has 
been explored and it has been found that negative dif- 
ferential conductance can be obtained for low values of 
the gate and bias applied voltages. Our results suggest 
that possible applications with GQDs can be developed 
for new electronic devices. 
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V. APPENDIX: GREEN'S FUNCTION OF THE 
SIMPLE MODEL 



In section III, we have introduced a simple model in 
order to explain the different slopes of figure |6] This 
model is composed by a linear chain of three sites, of the 
same energy, at which we have coupled four quantum 
dots (QDs) of energies £n {n = u,d), forming a crossbar 
junction configuration exhibited in figure [T] This simple 
scheme is very useful to qualitatively explain the elec- 
tronic behavior of graphene quantum dot that we have 
studied. 

Let us start with the hamiltonian of the system de- 
scribed by figure [7| 

Ht = Hieads Hc^leads (7) 

where the hamiltonian of the leads Hieads is given by: 



eads 



(8) 



the conductor hamiltonian He is given by; 



E E 

m( — 1,1) n{u,d) 



(9) 



and finally the leads-conductor hamiltonian is given by: 

Hc,leads= {fL Ck,a ^ he) (10) 

k,a{L,R) m(-l,l) 



By using the Dyson equation, it is possible to calcu- 
late the Green' function of site 0. Following a standard 
procedure we have obtained: 

Goo = go ^ govGol + govGol (11) 
Gio = givGoo ^ giYl ^R^k^^o (12) 

k 

Gio = govGio + ^i^VlG/c^,o, (13) 



where Gk^,o = gk VrGio and Gk^^o = gk VlGio- 

Replacing these expression in the previous set of equa- 
tions, we obtain: 



Gio = 
Gio = 



givG, 



00 



1 - gi Ekn 9kn 

gi V Goo 

1-^1 EkVlg,,VL 



(14) 
(15) 



Considering go = l/(cc; — ^o), and replacing the above 
expressions for Gio and Gjo, Goo reads: 



G, 



1 



00 



uj - eo 



(16) 



where the self-energy is defined by the following expres- 
sion: 



E = 



gi V 



gi ^ 



1-^1 ^r^^ 1-^1 i^L 



(17) 



with = E,, Vlgu, Vr and F^ = E,, vIqu^Vl. 

Using the expressions for the on-sites Green's functions 
for the sites 1 and 1 given by: gi = l/{uo — ei) and 
gi = l/(cj— £i), and considering a gate voltage Vg applied 
to the QDs, which redefine their on-sites energies by 

= ^n^ygi it is possiblc to Write an expression for the 
self-energy of the systems as: 



(18) 



where we have considered a symmetric system {Tl= Tr). 
In this approach, it is possible to write a compact form 
for the self-energy given in equation (18), which contains 



a real ( gate voltage dependent) and an imaginary part. 
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